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LaCu02.5 is a system of coupled, two-chain, cuprate ladders which may be doped systematically by 
Sr substitution. Motivated by the recent synthesis of single crystals, we investigate theoretically the 
possibility of superconductivity in this compound. We use a model of spin fluctuation-mediated su- 
perconductivity, where the pairing potential is strongly peaked at 7r in the ladder direction. We solve 
the coupled gap equations on the bonding and antibonding ladder bands to find superconducting 
solutions across the range of doping, and discuss their relevance to the real material. 



PACS numbers: 75.10.Jm, 75.40. Cx, 75.40. Gb 

In the search for the microscopic origin of high-T c su- 
perconductivity (HTS), it is essential to understand not 
only the materials in which it occurs, but also why it 
fails to occur in certain systems which are ostensibly 
only marginally different. In this letter we examine the 
cuprate Lai_ 2; Sr 2 ;Cu02.5, a coupled-chain analog of pla- 
nar La2-a;Sr a; Cu04. The parent compound was at first 
thought to be a dopable, two-chain spin ladder system 
, but it was later shown |^-^| that interladder coupling 
in this unfrustrated structure was sufficiently strong to 
stabilize long-range antiferromagnetic (AF) order. 

The first doping experiments ^]|| were on ceramic 
samples, and found metallic behavior only for relatively 
high doping values x > 0.2. No superconductivity was 
found for any doping, a result attributed to quasi-one- 
dimensionality (Id) and the unavoidably strong random 
potential ||. Recently, Takagi and coworkers have 
successfully grown single crystals with the fixed compo- 
sition x — 0.15. These show good metallic character- 
istics, both along and across the ladder direction with 
anisotropy ratio px/p\\ ~ 8, and residual resistivities as 
low as 40/if2cm~ 1 . Although the latter value corresponds 
to a mean free path exceeding 10 ladder unit cells, at least 
comparable to coherence lengths in cuprates, there was 
no indication of superconductivity. The subtlety of the 
HTS phenomenon was further highlighted by the discov- 
ery of superconductivity under pressure in a different 
cuprate ladder system. 

The resistivity data appear to support an 
anisotropic, 3d approach. Here we consider in detail the 
Fermi surface of Lai_ a: Sr a ;Cu02.5, to find out if there are 
inherent frustrations present which could hinder or pre- 
vent a superconducting state with the predominantly d- 
symmetric character expected in ladders || . To this end 
we use the model of Millis, Monien and Pines (MMP) 
Jill for superconductivity mediated by AF spin fluctu- 
ations, and perform exploratory calculations examining 
the form of the resulting gap function. 

Considering the structure as two-chain ladders along 



z, with weaker interladder couplings in the transverse di- 
rections, we compute tight-binding bands analogous to 
the bonding and antibonding combinations in an iso- 
lated ladder. In Ref. ||, a set of near-neighbor hop- 
ping parameters was deduced from a four-band fit to 
the band structure computed in the local density ap- 
proximation (LDA) ijyj. Using local, Cu-based orbitals 
4>i and 02 on each atom in the reduced unit cell, the 
Hamiltonian in terms of reciprocal-space wave functions 



6 M (k) = v- 1 ' 2 M^y k ' r ' (M = 1, 2) is 



til ti2 
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where 



tn = t 2 2 = -2t z cosk z 
ti2 = t 2 i = —t a — tb[e 



2t' z cos 2k z 
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U denotes t, + 2t[ cos k z , the ladder lattice constant a c 
has been set to 1, and the parameters ti, shown in Fig. 
1, have the values j| t z = t a = t = 0.4eV, t' z = \t Zl 
t' a = _ |iai t{, = ^t a and t' b = ^th- 




FIG. 1. Two-chain ladder structure of Cu ions in 
Lai_ a; Sr :E Cu02.5, showing tight-binding parameters. 

The two eigenmodes are £k± = in(k) ± |t 12 (k) | . The 
nature of the corresponding eigenfunctions is seen most 
clearly by introducing the bonding and antibonding basis 



1 



functions 0h iQ (k) = [0i(k)±</>2(k)]/\/2, in terms of which 
the eigenfunctions are 



: (k) = cos6> k Q:6 (k) + zsin0 k 0b,a(k), 



(4) 



with the phase factor #k =iarg(ti2). At the zone cen- 
ter these wave functions have full bonding or antibonding 
character, which is exchanged at the (tt, tt, 0) point. Else- 
where in the zone the character is mixed, except along 
the line k x = —k Vl where the argument remains zero and 
a band crossing occurs. The nature of the eigenfunctions 
changes little with k z , except close to (tt, tt, 0) where there 
is another crossing. The two-band dispersion for a dop- 
ing of x — 0.15 is illustrated in Fig. 2(a) for a k-space 
path equivalent to that depicted in Refs. |11| and ||. 
Although no crossings occur on this path, the character 
of the associated eigenfunctions is altered by evolution 
of the complex phase factor. The Fermi surfaces (Fig. 
2(b)) for the two bands remain quasi-ld ||, in that they 
take the form of sheets crossed only along k z , with some 
variation of k z -p depending on k = {k x , k y ). 
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FIG. 2. (a) Tight-binding dispersion relations and (b) 
Fermi surfaces for two-band model at doping x = 0.15. 

The gap equations describing superconductivity in a 
two-band system take the form 



A kQ = 77 > , Vkic'a/3— — tanh 



V ^ 

k',/3 



2E] 



k'/3 



2k B T 



(5) 



where a,/3 = ± are the band indices, Vkk'a/3 is a matrix 
containing both the on-site repulsion term and the pair- 
ing potential, and Ek a — ±y/(ek a — n) 2 + A kQ specifies 



the dispersion relations in the superconducting state (// 
is the chemical potential). The symmetry of the gaps 
Ak±, both along (k z ) and transverse (k) to the ladder 
axis, is of particular interest. 

The pairing potential, as usual in cuprates, is repulsive. 
Scattering processes due to interactions with spin fluctu- 
ations should be strongly enhanced at the characteristic 
wave vector of the AF correlated ladder spins, q z = n. 
The dominant contributions are expected from interband 
scattering of intraband electron pairs when the band sep- 
aration is close to 7r, as represented schematically in Fig 
3. To describe this we adopt the frequency-dependent 
MMP potential in the form 



Vkk'a/3 



-Vo6 a /3 



^(\+e z \q z -Q z \ 2 f + {^/^h) 



(6) 



Here q z = k z — k' z , Q z — tt is the AF wave vector in 
the ladder direction, ui = — e^/ and wsf sets the fre- 
quency scale. e Q( 3 is the off-diagonal tensor ensuring the 
interband scattering condition, and has the effect of a 
projector, Y, mi tn I0m)(0n|, m,n = a, 6, restricting the 
processes of interest to those scattering pairs between 
bonding and antibonding states. The longitudinal corre- 
lation length £ 2 determines the spin fluctuation enhance- 
ment, and may be deduced from the MMP parameters 
|]l0| , ^2| for planar cuprate superconductors. We neglect 
the transverse dispersion of spin fluctuations: because 
[Q,H interladder superexchange is weak (J± ~ 0.12 J), 
this would contribute an effect of order 10%, with the 
consequence that the leading components of any gap so- 
lution will be isotropic in k. 
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. 3. Schematic Fermi-surface geometry showin 

tant interband pair scattering processes. 
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The Hubbard-type model under consideration has also 
a repulsive on-site interaction U , which appears in the 
Hamiltonian as Htj =^Z)ka=±(^kQC k ]c"^ + H.c). 
In the appropriate (t-J) limit of U — > oo, the amplitude 
for pair creation on any site must be zero, which sets the 



constraint c^jc^j^ 



c k | c. 
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= 0. 



The full, coupled gap equations can be written as 



where 



2 Dk ' a 



Aka / i^k 

Oka = tt-^ — tanh 



2-Eka 



(7) 
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a = ±, a = —a, and i? denotes the denominator in 
Eq. (||). As both U and V terms are repulsive, the 
most straightforward solution is one where the gaps have 
opposite signs on the ± bands. The symmetry of the U 
terms in D + and D_ ensures that the constraint ^ is 
obeyed. 

We take the parameter values from planar cuprates 
fl2| ], where the local physics can be expected to be sim- 
ilar. The correlation length is £ z = 2.3o c . The prefactor 
is Vo = g 2 XQi where g = 1.53eV is the coefficient of the 
electron-spin fluctuation interaction and xq — 7T (^/ a ) Xo 
is a wave vector-adjusted static susceptibility; xo — 2-62 
states/eV gives V ~ 100. The damping frequency is 
given by ivsf = r/^C/a) 2 , where T = 0.4eV is the fre- 
quency cutoff of the interaction, whence usf = 8meV. 
The values of all parameters are taken to be constant at 
and below T c . Considering the range of U, the constraint 
(0) is satisfied to within 1% for U = 10 J, so this value is 
used below. 

Solution of the gap equations (|^) is simplified by sepa- 
rating the dependences of the pairing potential on k z and 
k, and a harmonic decomposition of the gap functions 



A 1 cos/cz 



A* cos 2k z 



and the potential 



V(k, k' z ) = Vy (k) cosik z cos jk' z . 



ij 



(10) 



(11) 



The dominant contributions to Vy(k) are from Fermi- 
surface scattering processes with 6k = 0. The approxi- 
mate treatment is justified by the parameters £ 2 and u>sf. 
which are such that the halfwidth of V(q z ) is Sq ~ 0.27T, 
a quantity exceeding the variation in band separation 
fei, — k~ F across the Brillouin zone (Fig. 2(b)). Thus 
scattering vectors of all q contribute similarly, and the 
averaging performed in the approximation introduces in- 
significant errors. 

The qualitative nature of the solutions is influenced by 
the doping. In this treatment the band shape is not af- 
fected, but a change of doping alters the separation of the 
Fermi surfaces (Fig. 3). The pairing is (broadly) peaked 
where the bands have maximal areas separated by q z ~ tt. 



We find that at x m = 0.112 the separation is 7r at k = 
(0,0), and is smaller than 7r everywhere else in the zone. 
Increased doping [x > x m ) causes the net interaction to 
diminish, while at x < x m a circle of band separation 
7r will open and increase in circumference, implying a 
monotonic rise in any T c , at least until x — 0.082 where 
Fermi surface is lost to formation of an electron pocket 
at Z in the — band. We have solved the gap equations 
for a variety of doping values 0.05 < x < 0.2, retaining 
interaction coefficients Vij up to fourth order, and gap 
components A, to i = 2. 




FIG. 4. Superconducting gaps (a) Ajt" and (b) on the 
Fermi surfaces for doping x = 0.082. 

A superconducting solution with finite gap components 
is obtained for all the doping levels considered. The gaps 
have opposite signs, and no nodes. Figs. 4(a) and (b) 
show these full gaps on the Fermi surfaces for doping 
x = 0.082, which is qualitatively representative of all 
cases. A^ is effectively constant, while has only ap- 
preciable zeroth and first harmonics, justifying the ne- 
glect of higher components. That the order parameters 
have opposite signs on each band corresponds within the 
ladder to the condition |ll| for "d-like" symmetry that 
A(fcj^ = 0) and A(fej_ = tt) have opposite signs. The 
order parameters are isotropic transverse to the ladder 
direction. The decrease of gap components with doping 
is quantified by T c (Fig. 5). The slower increase with de- 
creasing doping below x — 0.082 may be accounted for as 
above. Quantitatively, the values of T c predicted from the 
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two-band model with spin fluctuation-mediated pairing 
approach those of the 2d planar cuprate La2-a;Sr x Cu04. 
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FIG. 5. T c as a function of doping for two-band model. 

Our primary result is then that, within a specific 
model for AF spin fluctuations, the coupled ladder com- 
pound Lai_ x Sra;Cu02.5 is predicted to be superconduct- 
ing. Since the solutions of the gap equations have a full 
gap over the whole Fermi surface with opposite signs on 
the different sheets, there is clearly no inherent frustra- 
tion for a superconductor with approximate d-wave sym- 
metry. From this point of view, it appears rather that the 
Fermi surface form would favor a superconducting state. 
The primary weakness of the preceeding calculation is 
that it is carried out for a clean system, whereas in real- 
ity there are two sources of disorder. One is crystalline 
imperfections, which presumably may be controlled, and 
the other is intrinsic to the doping process which substi- 
tutes divalent Sr 2+ ions for trivalent La 3+ in close prox- 
imity to the ladders. Evidence for the disruptive influ- 
ence of this random potential may be found in recent ex- 
periments on thin films of Lai.gSro.iCu04 14|, in which 
the distance of Sr 2+ from the cuprate planes was varied. 
For greater separation, the transition temperature was 
substantially enhanced, from 29K to T c = 50K, and si- 
multaneously the in-plane resistivity in the normal state 
was reduced. This source of disorder cannot be elimi- 
nated, and may lead to localization, particularly at lower 
x, which would produce a doping dependence of T c of 
the type manifest in the planar cuprates (c/. Fig. 5). 
The synthesis of crystalline samples with x = 0.15 and 
low residual resistivities JtJ is encouraging in the follow- 
ing regard. It provides hope that if high-quality single 
crystals could be synthesized for smaller x, there is the 
possibility of superconductivity, or at minimum of the ob- 
servation of spin-gap features similar to the underdoped 
cuprates. If in the latter case a crossover to an insulat- 
ing phase appeared at low temperatures, this would be 
prime facie evidence that indeed the presence of disorder 
was responsible for a suppression of superconductivity. 

A related question concerns the Landau Fermi liquid 
nature of the metallic state observed at x = 0.15. At first 
sight, the quasi-ld nature of Lai_ x Sr x Cu02.5 should ac- 



celerate the breakdown of Landau theory, yet it seems 
that the critical doping required to stabilize it is actu- 
ally lower than in the planar cuprates. One possible 
explanation for this unexpected discrepancy is the role 
of Umklapp scattering, which has been proposed |l5j as 
the cause of the breakdown of Landau behavior in 2d 
cuprate systems. In the present case, Umklapp scatter- 
ing between the two Fermi-surface sheets occurs only for 
smaller values of x < x m = 0.112. Again it would be 
most interesting to follow the evolution of the normal 
state in samples with lower doping. 

In conclusion, we have found that from the point of 
view of spin-fluctuation theory, the Lai_ 2; Sr 2 ;Cu02.5 sys- 
tem, with its clear ladder structure, should favor a node- 
less, d-wave superconducting state. However, this may 
be suppressed by the strong intrinsic disorder introduced 
by substituting Sr 2+ for La 3+ . The synthesis of crys- 
talline samples with low doping values would shed light 
on many key issues in the understanding of cuprates. 
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